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, Abstract. Considerable progress has been made in understanding the 

lO ' hydrodynamics, but only to a certain extent the magnetohydrodynamics, 

of shaping bipolar outflows forming protoplanetary nebulae (PPNs) and 
planetary nebulae (PNs). In particular, Blackman et al. (2001b, 2001a) 
^ . point out two problems related to the formation of PNs and PPNs, regard- 

00 ' ing the formation of multipolar structures and the origin of the nebulae. 

0^ ■ They propose a solution by giving a semi-quantitative physical model 

I which should be investigated by numerical simulations. 



Motivation 



O I The Shape Problem. All current models involving magnetic fields assume 

' that a radiation driven fast post-AGB (asymptotic giant branch) wind carries 

! a (frozen-in!) stellar toroidal magnetic field, travelling into a non-magnetized 

radiation driven slow AGB wind. However, it is not clear whether this kind of 
magnetic shaping can be applied to nonaxisymmetric multipolar structures that 
^ ' have been observed by the Hubble Space Telescope. The multipolar structures 

suggest that two winds are interplaying, emanating at the same time from the 
star and a surrounding (accretion) disc, a scenario which has not been considered 
in previous models. 

The Power Problem. All current models assume that the winds are 
purely radiation driven. There is observational evidence that the PN-shaping 
process begins already in the PPN phase, because many PPNs have highly col- 
limated fast bipolar outflows. However, at this early stage of evolution, a post- 
AGB star is too cool to produce a high speed (radiation driven!) fast wind with 
the observed large kinetic luminosity. Therefore, the origin of fast PPN winds 
cannot be radiation, as assumed in previous models. 

The Proposed Model. Blackman et al. (2001b, 2001a) propose that large 
scale dynamos in the post-AGB star as well as a surrounding PPN accretion 
disc might drive MHD stellar and disc winds, so that PPNs/PNs are formed by 
dynamo driven stellar and disc winds rather than radiation driven stellar winds. 

Our Work. We aim, in the future, to verify the proposed semi-quantitative 
physical model with global numerical MHD simulations. As a first step, we 
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Figure 1. Model A. Left: Black arrows: poloidal velocity vectors; 
white solid lines: poloidal magnetic field lines; black dashed line: fast 
magnetosonic surface; colours: temperature (specific enthalpy h is di- 
rectly proportional to temperature T, and log^Q h = (—2, —1,0, 1) cor- 
responds to T (3x 10^, 3x 10^, 3x 10^ 3x 10^) K ). Right: Colours: 
specific angular momentum i (normalized with the maximum in the 
disc); white solid line: Alfven surface; white dotted line: sonic sur- 
face; white dashed line: fast magnetosonic surface. Both: Averaged 
over times t ^ 897 days . . .t ^ 906 days; the diffusion time is about 
576 days. 

present here results of axisymmetric numerical simulations of three MHD mod- 
els of a system comprised of a central object, a surrounding accretion disc and 
a corona. All three models have mean-field dynamo action in the disc and do 
not impose any external magnetic field in the disc, which is a new approach 
in the context of outflows and accretion studies. Model A assumes an initially 
non-magnetized star, whereas Models B and C include a strong stellar dipo- 
lar magnetosphere and a stellar mean-field dynamo, respectively. Our models 
are currently applied to protostellar star-disc systems. However, since we use 
dimensionless variables, the models can be rescaled and applied to a range of 
different astrophysical objects. 

2. Models and Results 

We solve the continuity equation, the Navier-Stokes equation and the induction 
equation. Instead of solving the energy equation, we model a cool dense disc 
embedded in a hot rarefied corona by fixing entropy contrasts such that spe- 
cific entropy is smaller within the disc and the star, and larger in the corona. 
By solving the hydrostatic equilibrium for an initially non-rotating corona, the 
corona is pressure-supported, whereas the disc is mainly centrifugally supported 
and rotating slightly sub-Keplerian. 

Model A: Accretion Disc Dynamo. The results of Model A have 
been presented in more detail in von Rekowski et al. (2003). As one can see 
in Fig. n (Left), a clearly structured outflow develops with a slower, hotter 
and denser stellar wind, an inner disc wind (in a conical shell originating from 
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Figure 2. Model B. Left: Same as in Fig. ^ but for Model B, and 
at a time (t ~ 158 days) when the star is about to be magnetically 
reconnected to the disc. Black dashed line: Alfven surface. The stellar 
surface magnetic field strength is ~ 5 kG. Middle and Right: Colours: 
density (bright colours indicate high values; dark colours indicate low 
values); white solid lines: poloidal magnetic field lines; black arrows: 
azimuthally integrated mass flux density (except in disc where density 
is high and vectors would be too long). The stellar surface magnetic 
field strength is ~ 5 kG. Middle: At a time {t ~ 150.5 days) when 
the star is magnetically connected to the disc. Right: At a time (t ~ 
156.5 days) when the star is magnetically disconnected from the disc. 



the inner disc edge) that is faster, cooler and less dense, and an outer disc 
wind with intermediate values. The stellar wind is pressure driven (very low 
angular momentum, see Fig. ^ (Right); ratio of the poloidal magneto-centrifugal 
to pressure forces is less than 1 at the stellar surface) , whereas the inner disc wind 
is magneto-centrifugally accelerated (very high angular momentum; magnetic 
field lines at the disc surface are inclined to the rotation axis by at least 30°; 
outflow is highly supersonic but sub-Alfvenic: magnetic lever arm is about 3; 
ratio of the poloidal magneto-centrifugal to pressure forces is larger than 1 at 
the disc surface and exceeds 100 further away), and the outer disc wind is mostly 
pressure driven (Alfven surface is close to the disc surface; ratio of the poloidal 
magneto-centrifugal to pressure forces is mostly less than or around 1 at the disc 
surface). The averaged accretion rate is about 10~^Mo/yr and the averaged disc 
+ stellar wind mass loss rate is about 4 x lO^''M0/yr. 

Model B: Accretion Disc Dynamo & Stellar Magnetosphere. The 
results of Model B have been presented in more detail in von Rekowski & Bran- 
denburg (2003a). Again, the disc dynamo produces a structured disc wind so 
that the overall outflow is similar to Model A, with coexisting pressure driving 
and magneto-centrifugal acceleration. In transition periods when disconnected 
stellar and disc field lines are about to reconnect, there is a pressure driven 
hot and dense, but relatively fast outflow between the stellar and disc winds 
(see Fig. |21 Left). The magnetosphere is oscillating with a period of around 
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Figure 3. Model C. Same as in Figs EandEl (Left), but for Model C, 
and at a time t 1172 days. 

15 to 30 days, changing between two distinct states. When star and disc are 
magnetically connected by the magnetosphere, accretion of disc matter onto 
the star is along magnetospheric field lines with peak accretion rates of about 
(1 . . . 2.5) X lO~^M0/yr (see Fig. |2J Middle). When star and disc are magneti- 
cally disconnected, there is no net accretion but matter is lost directly into the 
outflow (see Fig. [2 Right). In this state, angular momentum transport by the 
inner disc wind is clearly enhanced; also the stellar and inner disc wind veloc- 
ities are enhanced. The disc wind mass loss rates are also time dependent but 
not correlated with the oscillating magnetosphere; they are about one order of 
magnitude higher than the peak accretion rates. 

Model C: Accretion Disc Dynamo & Stellar Dynamo. The results 
of Model C are going to be presented in more detail in von Rekowski & Branden- 
burg (2003b). The stellar dynamo generates a magnetic field with dipolar- type 
symmetry (the alpha effect is positive in the upper hemisphere). The main re- 
sult is that the dynamo in the star produces a fast and collimated stellar wind 
(cf. Fig. 131). The structured disc wind is again due to the disc dynamo. 
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